Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Ito, Takashi; Kadono, Ryosuke*
Journal of the Physical Society of Japan, 93(4), p.044602_1 - 044602_7, 2024/04
Riyana, E. S.; Okumura, Keisuke; Sakamoto, Masahiro; Matsumura, Taichi; Terashima, Kenichi; Kanno, Ikuo
Journal of Nuclear Science and Technology, 8 Pages, 2024/00
Times Cited Count:0 Percentile:0.18(Nuclear Science & Technology)Tuya, D.; Nagaya, Yasunobu
Journal of Nuclear Engineering (Internet), 4(4), p.691 - 710, 2023/11
The Monte Carlo method is used to accurately estimate various quantities such as k-eigenvalue and integral neutron flux. However, when a distribution of a quantity is desired, the Monte Carlo method does not typically provide continuous distribution. Recently, the functional expansion tally and kernel density estimation methods have been developed to provide continuous distribution. In this paper, we propose a method to estimate a continuous distribution of a quantity using artificial neural network (ANN) model with Monte Carlo-based training data. As a proof of concept, a continuous distribution of iterated fission probability (IFP) is estimated by ANN models in two systems. The IFP distributions by the ANN models were compared with the Monte Carlo-based data and the adjoint angular neutron fluxes by the PARTISN code. The comparisons showed varying degrees of agreement or discrepancy; however, it was observed that the ANN models learned the general trend of the IFP distributions.
Wakui, Takashi; Takagishi, Yoichi*; Futakawa, Masatoshi
Materials, 16(17), p.5830_1 - 5830_16, 2023/09
Times Cited Count:0 Percentile:0(Chemistry, Physical)Cavitation damage on the mercury target vessel is induced by proton beam injection in mercury. The prediction method of the cavitation damage using Monte Carlo simulations was proposed taking into account of the uncertainties of the position of cavitation bubbles and impact pressure distributions. The distribution of impact pressure attributed to individual cavitation bubble collapsing was assumed to be the Gaussian distribution, and the probability distribution of the maximum value of impact pressures was assumed to be three kinds of distributions; the delta function, the Gaussian and Weibull distributions. Two parameters were estimated using Bayesian optimization by comparing the distribution of the cavitation damage obtained from experiment with that of accumulated plastic strain obtained from the simulation. It was found that the results obtained using the Weibull distribution reproduced the actual cavitation erosion phenomenon better than the other results.
Kubo, Kotaro; Tanaka, Yoichi; Hakuta, Yuto*; Arake, Daisuke*; Uchiyama, Tomoaki*; Muramatsu, Ken
Mechanical Engineering Journal (Internet), 10(4), p.23-00051_1 - 23-00051_17, 2023/08
The significance of probabilistic risk assessments (PRAs) of nuclear power plants against external events was re-recognized after the Fukushima Daiichi Nuclear Power Plant accident. Regarding the seismic PRA, handling correlated failures of systems, components, and structures (SSCs) is very important because this type of failure negatively affects the redundancy of accident mitigation systems. The Japan Atomic Energy Research Institute initially developed a fault tree quantification methodology named the direct quantification of fault tree using Monte Carlo simulation (DQFM) to handle SSCs' correlated failures in detail and realistically. This methodology allows quantifying the top event occurrence probability by considering correlated uncertainties related to seismic responses and capacities with Monte Carlo sampling. The usefulness of DQFM has already been demonstrated. However, improving its computational efficiency would allow risk analysts to perform several analyses. Therefore, we applied quasi-Monte Carlo and importance sampling to the DQFM calculation of simplified seismic PRA and examined their effects. Specifically, the conditional core damage probability of a hypothetical pressurized water reactor was analyzed with some assumptions. Applying the quasi-Monte Carlo sampling accelerates the convergence of results at intermediate and high ground motion levels by an order of magnitude over Monte Carlo sampling. The application of importance sampling allows us to obtain a statistically significant result at a low ground motion level, which cannot be obtained through Monte Carlo and quasi-Monte Carlo sampling. These results indicate that these applications provide a notable acceleration of computation and raise the potential for the practical use of DQFM in risk-informed decision-making.
Li, C.-Y.; Wang, K.*; Uchibori, Akihiro; Okano, Yasushi; Pellegrini, M.*; Erkan, N.*; Takata, Takashi*; Okamoto, Koji*
Applied Sciences (Internet), 13(13), p.7705_1 - 7705_29, 2023/07
Times Cited Count:1 Percentile:0(Chemistry, Multidisciplinary)Wakui, Takashi; Takagishi, Yoichi*; Futakawa, Masatoshi; Tanabe, Makoto*
Jikken Rikigaku, 23(2), p.168 - 174, 2023/06
Cavitation damage on the inner surface of the mercury target for the spallation neutron source occurs by proton bombarding in mercury. The prediction method of the cavitation damage using Monte Carlo simulations was suggested taking variability of the bubble core position and impact pressure distribution into account. The impact pressure distribution was estimated using the inverse analysis with Bayesian optimization was conducted with comparison between cavitation damage distribution obtained from experiment and the cumulative plastic strain distribution obtained from simulation. The average value and spread of maximum impact pressure estimated assuming the Gaussian distribution were 3.1 GPa and 1.2 m, respectively. Simulation results reproduced experimental results and it can be said that this evaluation method is useful.
Li, C.-Y.; Watanabe, Akira*; Uchibori, Akihiro; Okano, Yasushi
Proceedings of 30th International Conference on Nuclear Engineering (ICONE30) (Internet), 10 Pages, 2023/05
Rodriguez, D.; Rossi, F.
Proceedings of INMM & ESARDA Joint Annual Meeting 2023 (Internet), 9 Pages, 2023/05
Kubo, Kotaro; Jang, S.*; Takata, Takashi*; Yamaguchi, Akira*
Journal of Nuclear Science and Technology, 60(4), p.359 - 373, 2023/04
Times Cited Count:5 Percentile:78.52(Nuclear Science & Technology)Probabilistic risk assessment (PRA) is an essential approach to improving the safety of nuclear power plants. However, this method includes certain difficulties, such as modeling of combinations of multiple hazards. Seismic-induced flooding scenario includes several core damage sequences, i.e., core damage caused by earthquake, flooding, and combination of earthquake and flooding. The flooding fragility is time-dependent as the flooding water propagates from the water source such as a tank to compartments. Therefore, dynamic PRA should be used to perform a realistic risk analysis and quantification. This study analyzed the risk of seismic-induced flooding events by coupling seismic, flooding, and thermal-hydraulics simulations, considering the dependency between multiple hazards explicitly. For requirements of safety improvement, especially in light of the Fukushima Daiichi Nuclear Power Plant accident, sensitivity analysis was performed on the seismic capacity of systems, and the effectiveness of alternative steam generator injection by a portable pump was estimated. We demonstrate the use of this simulation-based dynamic PRA methodology to evaluate the risk induced by a combination of hazards.
Li, C.-Y.; Watanabe, Akira*; Uchibori, Akihiro; Okano, Yasushi
Journal of Nuclear Science and Technology, 23 Pages, 2023/00
Times Cited Count:1 Percentile:72.91(Nuclear Science & Technology)Tuya, D.; Nagaya, Yasunobu
Nuclear Science and Engineering, 15 Pages, 2023/00
Times Cited Count:0 Percentile:0.01(Nuclear Science & Technology)In Monte Carlo neutron transport calculations for local response or deep penetration problems, some estimation of an importance function is generally required in order to improve their efficiency. In this work, a new recursive Monte Carlo (RMC) method, which is partly based on the original RMC method, for estimating an importance function for local variance reduction (i.e., source-detector type) problems has been developed. The new RMC method has been applied to two sample problems of varying degrees of neutron penetrations, namely a one-dimensional iron slab problem and a three-dimensional concrete-air problem. The biased Monte Carlo calculations with variance reduction parameters based on the obtained importance functions by the new RMC method have been performed to estimate detector responses in these problems. The obtained results are in agreement with those by the reference unbiased Monte Carlo calculations. Furthermore, the biased calculations offered an increase in efficiency on the order of 1 to 10 in terms of the figure of merit (FOM). The results also indicated that the efficiency increased as the neutron penetration became deeper.
Simanullang, I. L.*; Nakagawa, Naoki*; Ho, H. Q.; Nagasumi, Satoru; Ishitsuka, Etsuo; Iigaki, Kazuhiko; Fujimoto, Nozomu*
Annals of Nuclear Energy, 177, p.109314_1 - 109314_8, 2022/11
Times Cited Count:0 Percentile:0.01(Nuclear Science & Technology)Omer, M.; Shizuma, Toshiyuki*; Hajima, Ryoichi*; Koizumi, Mitsuo
Dai-43-Kai Nihon Kaku Busshitsu Kanri Gakkai Nenji Taikai Kaigi Rombunshu (Internet), 3 Pages, 2022/11
Zheng, X.; Tamaki, Hitoshi; Takahara, Shogo; Sugiyama, Tomoyuki; Maruyama, Yu
Proceedings of Probabilistic Safety Assessment and Management (PSAM16) (Internet), 10 Pages, 2022/09
Omer, M.; Shizuma, Toshiyuki*; Hajima, Ryoichi*; Koizumi, Mitsuo
Radiation Physics and Chemistry, 198, p.110241_1 - 110241_7, 2022/09
Times Cited Count:2 Percentile:53.91(Chemistry, Physical)Kubo, Kotaro; Fujiwara, Keita*; Tanaka, Yoichi; Hakuta, Yuto*; Arake, Daisuke*; Uchiyama, Tomoaki*; Muramatsu, Ken*
Proceedings of 29th International Conference on Nuclear Engineering (ICONE 29) (Internet), 8 Pages, 2022/08
After the Fukushima Daiichi Nuclear Power Plant accident, the importance of conducting probabilistic risk assessments (PRAs) of external events, especially seismic activities and tsunamis, was recognized. The Japan Atomic Energy Agency has been developing a computational methodology for seismic PRA, called the direct quantification of fault tree using Monte Carlo simulation (DQFM). When appropriate correlation matrices are available for seismic responses and capacities of components, the DQFM makes it possible to consider the effect of correlated failures of components connected through AND and/or OR gates in fault trees, which is practically difficult when methods using analytical solutions or multidimensional numerical integrations are used to obtain minimal cut set probabilities. The usefulness of DQFM has already been demonstrated. Nevertheless, a reduction of the computational time of DQFM would allow the large number of analyses required in PRAs conducted by regulators and/or operators. We; therefore, performed scoping calculations using three different approaches, namely quasi-Monte Carlo sampling, importance sampling, and parallel computing, to improve calculation efficiency. Quasi-Monte Carlo sampling, importance sampling, and parallel computing were applied when calculating the conditional core damage probability of a simplified PRA model of a pressurized water reactor, using the DQFM method. The results indicated that the quasi-Monte Carlo sampling works well at assumed medium and high ground motion levels, importance sampling is suitable for assumed low ground motion level, and that parallel computing enables practical uncertainty and importance analysis. The combined implementation of these improvements in a PRA code is expected to provide a significant acceleration of computation and offers the prospect of practical use of DQFM in risk-informed decision-making.
Li, C.-Y.; Watanabe, Akira*; Uchibori, Akihiro; Okano, Yasushi
Dai-26-Kai Doryoku, Enerugi Gijutsu Shimpojiumu Koen Rombunshu (Internet), 4 Pages, 2022/07
Identifying accident scenarios that could lead to severe accidents and evaluating their frequency of occurrence are essential issues. This study aims to establish the methodology of the dynamic Probabilistic Risk Assessment (PRA) for sodium-cooled fast reactors that can consider the time dependency and the interdependence of each event. Specifically, the Continuous Markov chain Monte Carlo (CMMC) method is newly applied to the SPECTRA code, which analyzes the severe accident conditions of nuclear reactors, to develop an evaluation methodology for typical external hazards. Currently, a fault-tree model of air coolers of decay heat removal system is implemented as the CMMC method, and a series of preliminary analysis of the plant's transient characteristics under the scenario of volcanic ashfall has been conducted.
Tuya, D.; Nagaya, Yasunobu
Annals of Nuclear Energy, 169, p.108919_1 - 108919_9, 2022/05
Times Cited Count:1 Percentile:31.61(Nuclear Science & Technology)Estimating an effect of a perturbation in a fissile system on its -eigenvalue requires special technique called perturbation theory when the considered perturbation is small. In this study, we develop an adjoint-weighted correlated sampling (AWCS) method based on the exact perturbation theory without any approximation by combining the correlated sampling (CS) method with iterated-fission probability (IFP) based adjoint-weighting method. With the advantages of the CS method being good at providing very small uncertainty for small perturbations and the IFP-based adjoint-weighting method being suitable for continuous-energy Monte Carlo calculation, the developed AWCS method based on the exact perturbation theory offers a new rigorous approach for perturbation calculations. The obtained results by the developed AWCS method for verification problems involving Godiva and simplified STACY density perturbations showed good agreement with the reference calculations.
Riyana, E. S.; Okumura, Keisuke; Sakamoto, Masahiro; Matsumura, Taichi; Terashima, Kenichi
Journal of Nuclear Science and Technology, 59(4), p.424 - 430, 2022/04
Times Cited Count:1 Percentile:16.35(Nuclear Science & Technology)